Phase diagram of one-dimensional Hubbard-Holstein model at quarter-filling 
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We derive an effective Hamiltonian for the one-dimensional Hubbard-Holstein model, valid in a 
regime of both strong electron-electron (e-e) and electron-phonon (e-ph) interactions and in the non- 
adiabatic limit {t/iuo < 1), by using a non-perturbative approach. We obtain the phase diagram at 
quarter-filling by employing a modified Lanczos method and studying various density-density corre- 
lations. The spin-spin AF (antiferromagnetic) interactions and nearest-neighbor repulsion, resulting 
from the e-e and the e-ph interactions respectively, are the dominant terms (compared to hopping) 
and compete to determine the various correlated phases. As e-e interaction (U /t) is increased, the 
system transits from an AF cluster to a correlated singlet phase through a discontinuous transition 
at all strong e-ph couplings 2 < gr < 3 considered. At higher values of U /t and moderately strong 
e-ph interactions (2 < g < 2.6), the singlets break up to form an AF order and then to a paramag- 
netic order all in a single sublattice; whereas at larger values of g (> 2.6), the system jumps directly 
to the spin disordered charge-density-wave (CDW) phase. 

PACS numbers: 71.10.Fd, 71.38.-k, 71.45.Lr, 74.25.Kc 
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I. INTRODUCTION: 

A host of materials show evidence of e-ph interactions 
besides the expected e-e interactions. Angle-resolved 
photocmission spectroscopy (ARPES) experiments in 
cupratcsii^. fuUcridcs^, and manganitcsA indicate strong 
e-ph coupling. The interplay of c-c and c-ph interac- 
tions in these correlated systems leads to coexistence of 
or competition between various phases such as supercon- 
ductivity, CDW, spin-density-wave (SDW) phases, or for- 
mation of novel non-Fermi liquid phases, polarons, bipo- 
larons, etc. It is of particular interest to consider both 
the strong coupling regime, and the cross-over to the non- 
adiabatic limit where the lattice response is not slower in 
comparison to the heavy effective mass of the correlated 
electron system. 

The simplest framework to analyze the effects of e-e 
and the concomitant c-ph interactions is offered by the 
Hubbard-Holstein model whose Hamiltonian is given by: 



hh = 



JO" 



'^]+lcr'^io" 



H.c. 



-Wo 



Uj^^jt^ji- (1) 



Here, 



'JO 



-jcr'-jo 



with Cj^ and a] being the creation op- 
erators at site j for an electron with spin a and a phonon 
respectively. The Hamiltonian describes a tight-binding 
model with hopping amplitude t, a set of independent 
oscillators characterized by a dispersionless phonon fre- 
quency ujQ, along with an onsite Coulomb repulsion of 
strength U and an onsite electron-phonon interaction of 
strength g. The Holstein and the Hubbard models are 
recovered in the limits [/ = and 5 = respectively. 

To gain insight into the rich physics of the Hubbard- 
Holstein model, several studies have been conducted (in 
one-, two-, and infinite-dimensions and at various fill- 
ings) by employing various approaches such as quantum 



Monte Carlo (QMC)^"— , exact diagonalizationii"— , den- 
sity matrix renormalization group (DMRG)^, dynami- 
cal mean field theory (DMFT)^"— , semi-analytical slave 
boson approximations^^"— , large-N expansion^^, varia- 
tional methods based on Lang-Firsov transformational, 
and Gutzwiller approximation^. 

In contrast to earlier approaches, we utilize a controlled 
analytic approach (that takes into account quantum 
phonons). Our method uses both the strong electron- 
phonon coupling limit g > 1 and the strong Coulomb 
coupling limit U/t > 1 to generate an effective t ~ J 
model with displaced oscillators that can then be treated 
perturbatively, and generates longer-range interactions 
in an effective Hamiltonian. This model we then solve 
numerically for finite chains. We then obtain the phase 
diagram of the Hubbard-Holstein model at quarter-filling 
in one dimension. 

Our effective Hamiltonian comprises of two dominant 
competing interactions - spin-spin AF interaction and 
nearest-neighbor (NN) electron repulsion. In addition, 
three types of hopping also result ~ the NN hopping with 
reduced band width, next-nearest-neighbor (NNN) hop- 
ping, and NN spin-pair aa hopping. As the e-e interac- 
tion U/t is increased, the system sequentially transforms 
from an AF cluster phase to a correlated singlet phase 
followed by CDW phase(s) with (e-ph coupling g depen- 
dent) accompanying spin order. The most interesting 
feature is that, at intermediate values of U/t and for all 
strong e-ph couplings (2 < g < 3) considered, a phase 
comprising of correlated NN singlets appears which sug- 
gests the possibility for superconductivity occurrence. 

The paper is organized as fallows. In section II, we de- 
scribe our non-perturbative approach for tackling strong 
e-ph interaction and derive the effective Hamiltonian; in 
section HI, we discuss the region of parameter space for 
U/t, LOo/t, and g where our derived results are applica- 
ble; next, different correlation functions are analyzed to 
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obtain the various phases and the conditions for transi- 
tions between these phases are presented in section IV; 
then, in section V, we present the phase diagram; and in 
section VI our conclusions. 



II. EFFECTIVE HAMILTONIAN 

To get the effective Hubbard-Holstein Hamilto- 
nian, we first carry out the weh-known Lang-Firsov 



(LF) transformatioi 



ttLF 

^hh 



e HhhG where 



T — ^g^j^njai^CLj — a])- This transformation clothes 
the hopping electrons with phonons and displaces the 
simple harmonic oscillators. On rearranging the various 
terms, we get the following LF transformed Hamiltonian: 



JO" 



^ n, + {U - 2.g2a;o) ^ '^jf'^ji, (2) 



where Xj = e^'°^ "j-* and Uj = njA;+nj^. On noting that 
XjXj = 1, we can express, as shown below, our LF trans- 
formed Hamiltonian in terms of the composite fermionic 
operator d^- = cJg.X,J (i.e., a fermionic operator dressed 



with phonons): 



+ (C/ - 2.g2c.o) ^^'Ai - ff'^o E ("^T + O ' (3) 

where rt^^. = dj^dja- The last term in Eq. ([3]) represents 
the polaronic energy and is a constant for a given number 
of particles and an associated set of parameters. Hence, 
we drop this term from now onwards. 

The model represented by Eq. ^ is essentially the 
Hubbard Model for composite fermions with Hubbard 
interaction Ueff = {U — 2g^u!Q) and can be directly con- 
verted to an effective t — J model Hamiltonian in the limit 
of large Ueff/t. Thus, we get the effective t — J Hamil- 
tonian for composite fermions by projecting out double 
occupation. 



Ht-j = Ps 



JO" 



-J«J+1 



(4) 



wher 



J = 



4t^ 
U-2g^LOo 



, Si is the quantum 



"3 ~ "jt ^ "ji' 
mechanical spin operator for a spin 1/2 fermion at site 
I, and Ps is the operator that projects onto the singly 



occupied subspace. Next, we note that 

— g-9^g9(a]+i-a])g-s(aj + i-Oj)^ ^g^j 

The effective t — J Hamiltonian, given in Eq. ^ , can 
be re-expressed in terms of fermionic operators as 



Ht-.i — Ho + Hi, 



where 

Ho = -te-s' 



JO 

J^Ps (Sj ■ Sj+1 — ^-j-^ 



Ps 



(6) 



(7) 



and 
Hi 



-te-3 



Here, we have separated the Hamiltonian into (i) an 
itinerant electronic system represented by Hq contain- 
ing nearest-neighbor hopping with a reduced ampli- 
tude (te~^ ), electronic interactions, and no electron- 
phonon interaction; and (ii) the remaining part Hi 
which is a perturbation and corresponds to the compos- 
ite fermion terms containing the e-ph interaction with 

A. Perturbation Theory 

The unperturbed Hamiltonian Hq is characterized by 
the eigenstates \n, m) = \n)ei ® \'m)ph and corresponding 
eigencnergies En^m = En'^P'ru ■ noting that the first- 
order perturbation term is zero (i.e., (0, 0|i/i |0, 0) = 0), 
we proceed to calculate the second-order pertur- 
bation term E^^) = V (o.o|giK»»>(».rn|gi|o,o) ^ 

n,rn 

manner similar to that introduced in Ref. [sol Now, 
Ai?„j = E^l^ — E^^ is a positive integral multiple of 
ujQ and £"„ — Eq ^ te~3 = teff. In the narrow band 
limit, that te// <^ ujq <^ Heff, we get the corresponding 
second-order perturbation term in the effective Hamilto- 
nian for the polarons to be (see Appendix A for details) : 



ph{0\ Hi \m)ph X ph{m | Hi 1 0)p/t 



(9) 



Evaluation of ij'-^^ leads to the following expression: 

^ Y Ps [V {cl^i,c,,cl,c^^iy + c^c^+iJ^^i^^c^.' 



JOCT 

+ t2 {c]_. 



Ps, (10) 
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where V = (2/i + f2)tlff/ujo and t 

OO r, OO . 

A = E £r and E ^ 



2 — 

2(m + n) 

n!m! (m+n) 



fitljf/ujQ with 
. The value of 

/t;//t;i//t-t-/ti 

n.rn— 1 

/i and 2/i + /2 can be approximated for large value of 
g as / and e^^ /2(7^ respectively. But in numerical 
simulation we have calculated the value of fi and /2 by 
summing over the actual series. 

Finally, to calculate each term in Eq. pO)). we must 
project out the double occupancy. We do this projec- 



tion by replacing every fermionic operator Cj^ with the 
fermionic operator Cja-{1 — rijg) (see Appendix B for de- 
tails). 



B. Effective Electronic Hamiltonian 

The effective Hamiltonian, after averaging over the 
phononic degrees of freedom, is given by 



^hh — 



^(1 ~nj+is){l -njs){nja - rij+i^f 



^(1 - nj+is){l - nja){l - rij^ia) c]^i^{l - 2nj^)cj-ia + H.c 



JO" 



(11) 



JO" 



Here, we have approximated the coefficients as ^ ~ 
t'^/2g'^uJo and t2 — t'^e~^' / g'^uJo which is valid for large 
values of g. The operator form of the last three terms 
in Eq. (jlip can be visualized by considering different 
hopping processes depicted in Fig. [T] The third term 
in Eq. (fTTj) depicts the process where an electron hops 
to its neighboring site [see Figs. [T{a) andlljb)] and re- 
turns back. These two processes add up to the term 
(rijo- — nj+icr)^ whose expansion yields the NN repulsion 
term 2nj^nj^i^. The fourth term is composed of the two 
hopping processes shown in Figs. [TJc) andjUd). Fig. [TJc) 
represents a double hopping process where a spin a par- 
ticle hops first to its NN site and then to its NNN site. 
Contrastingly, Fig. [T][d) depicts the sequential process 
involving a pair of spin a electrons where an electron at 
site j hops to its neighboring site j + 1 followed by another 
electron at site j — I hopping to site j. This sequential 
process may be called era pair hopping. The fifth term 
corresponds to Fig. [T](e) and represents a hopping process 
similar to that shown in Fig. [Ijd) but involving a pair of 
electrons with opposite spins acr. The additional factors 
(l-nj+iff)(l-njff), (l-nj+iff)(l-njff)(l-nj_i5), and 
(1 — — nj-icr) appearing respectively in the third, 
fourth, and fifth terms in Eq. (|f ip represent projecting 
out double-occupancy at a site. The detailed derivation 
of these projected terms is given in Appendix [Bl 

The coefficients of the last three terms in Eq. (fTT|) (i.e., 
t^/2g^a;o and t^e~^ /g'^ujo), obtained from second or- 
der perturbation theory, are explained with the help of 
schematic diagrams shown in Fig. [2] We consider two 
distinct time scales for electronic hopping processes be- 
tween two adjacent sites: (i) ~ t~}f associated with either 



full distortion of the lattice ions to form a small pola- 
ronic potential well (of energy —g'^ujo) or full relaxation 
from the small polaronic distortion and (ii) ~ dur- 
ing which negligible distortion/relaxation of lattice ions 
occurs. The upper process in Fig. O^a) [schematically 
representing the distortion for hopping sequence in Fig. 
[Ija)] shows the process in which an electron goes to its 
right neighbouring site and comes back. The interme- 
diate state has energy +g'^u!Q and corresponds to fully 
distorted site j but without the electron and neighbour- 
ing site j + I containing the electron but without lattice 
distortion. To go from the intermediate state to the final 
one, the electron at site j + 1 hops back to the original 
site j without any new distortion taking place. Thus, 
the initial, the intermediate, and the final states all have 
identical lattice distortions. Hence, the hopping times 
for both these processes is t~^. Then, from second or- 
der perturbation theory, the numerator of the coefficient 
is while the denominator (which is energy difference 
between the intermediate and the initial states) becomes 
2(7^wo leading to the coefficient jpZJ^- ^^'^ same coeffi- 
cient results from the process in which an electron goes 
to its left neighbouring site and comes back as depicted 
by the lower process in Fig. [SJa). Thus the upper and 
lower processes in Fig. ^a.) both yield the same coeffi- 



The coefficient of the fourth term t^e~^ /s^^o corre- 
sponds to the schematic distortion processes shown in 
Figs. [2Ib) andlljc) with the pertinent hopping processes 
being depicted in Figs. [TJc) and [TJd) respectively. For 
the process where one electron at site j — 1 consecutively 
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hops to its NN site j and then to NNN site j + 1, the in- 
termediate states which give dominant contributions are 
shown in Fig. EJb). After the electron hops from site 
j — 1 to site j, the upper intermediate state in Fig. [2{b) 
has the same lattice distortion as the initial state. On the 
other hand, again in the upper intermediate state in Fig. 
m^b), when the electron makes the next hop from site j 
to _7 + 1 to produce the final state, there is a distortion 
at site J + 1 with a concomitant relaxation at the initial 
site j — I. Hence the coefficient contribution, from the 
hopping process involving the upper intermediate state 
in Fig. EJ^b), becomes t x te~^ /2g'^uJo. Alternately, af- 
ter the electron hops from site j — 1 to site j, the initial 
state may lead to the lower intermediate state in Fig. 
[2ljb) where the distortion at site j — 1 is completely re- 
laxed while site j + 1 gets distorted simultaneously so 
that the intermediate and the final states have identical 
lattice distortions. This process too yields a coefficient 
contribution of te~^ x t/2g^ujo. The other six possible 
intermediate states (corresponding to the remaining pos- 
sibilities of full distortion/relaxation at the three sites) 
each give t^e~^^ in the numerator and hence we ig- 
nore them. 

The coefficient of aa or ctct pair hopping [depicted by 
Figs. [Ijd) and [Ije) respectively] can be derived from 
Fig. [2Jc) which schematically represents the dominant 
contributions. The hopping process in Fig. [2l^c), with 
upper (lower) intermediate state, represents sequential 
hopping where an electron at site j hops to + 1 and 
produces the following changes: site j -I- f is undistorted 
(distorted), j — 1 distortion is unchanged (relaxed), and 
j remains distorted. The hopping time for the first hop 
is \/t (l/teff) and the change in energy between inter- 
mediate and initial states is 2g^u!o. Next, the electron at 
site J — 1 hops to j with site j — 1 relaxing (remaining 
undistorted) and site j + 1 getting (remaining) distorted. 
The second hop occurs in time 1/te// (I/O- Thus the aa 
or a<7 pair hopping process, involving the two dominant 
intermediate processes of Fig. ^c), yields the coefficient 



III. REGION OF VALIDITY OF OUR THEORY 

Our theory has three dimensionless parameters - Hub- 
bard interaction U/t, adiabaticity t/ioo, and electron- 
phonon interaction strength g. An increase in U/t, while 
keeping fixed other dimensionless parameters, will de- 
crease only the spin Hcisenberg AF interaction J = 
u-2g'^ujo ^'^ effective Hamiltonian of Eq. Fur- 
thermore, an increase in g produces an increase in J, 
a decrease in V, and also a decrease in all the various 
hopping (i.e., NN hopping, double hopping, and pair 
hopping) coefficients in Eq. ([TT]) . 

For large values of g, the coefficients of NN hopping, 
NNN hopping, {era) pair hopping terms contain and 
hence are small compared to NN repulsion V and spin in- 








j-1 j j+1 

(e) 

FIG. 1. Different hopping processes contributing to sec- 
ond order perturbation theory are: (a) c^lj^Cj+iCTcj^j^^Cjo-; (b) 
'-j+io-'-iCTcj^Cj+icr ; (c) cj^jjjCjo-cj^Cj_iCT ; (d) cj^Cj_iCTcj^-^^CjCT; 
and (e) c^j^Cj^xaC^jj^-^c^Cja. Here, empty circles correspond to 
sites with no electrons while circles with big and small dots 
inside correspond to sites with spin g and spin g electrons 
respectively. The numbers 1 and 2 indicate the order of hop- 
ping. 



teraction J. Thus, unlike the usual t—J model, here the 
spin-spin interaction term dominates over the hopping. 
Our model basically describes the competition between 
spin interactions and NN repulsion with the hopping be- 
ing a small perturbation. By varying the dimensionless 
parameters in our model, the relative importance of J 
and V terms can be changed. When the effect of J dom- 
inates over that oi V we expect formation of an AF 
cluster. On the other hand, when NN repulsion is larger 
than spin interaction, the electrons tend to get separated 
which can lead to different interesting phases. 

For our theory to be valid, the following criterion need 
to be satisfied: (a) the dimensionless effective Hubbard 
interaction Ueff/t ^ 1 so that Eq. ([3]) can be well ap- 
proximated by the effective t ~ J model given in Eq. 
(|4|). Actually, for a one-dimensional Hubbard model, 
calculations^ show that Ueff/t > 8.0 is sufficient to re- 
move double occupancy; (b) NN electronic hopping Ji 
and spin interaction J should both be negligible com- 
pared to the phononic energy wq so as to make our pertur- 
bation theory valid; and (c) the small parameter— of 
our perturbation theory t/ {gujo) should be kept as small 
as possible. Our calculations are done for t/ gu^ < 0.5. 
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(a) \®/ . \w7 '.'J 



=> 



j j+i 



j j+i 



j i+1 




j-i j j+i 

Initial state Relevant intermediate states Final state 

FIG. 2. Schematic representation of the hopping processes 
considered in Fig. [l]with intermediate states giving the dom- 
inant contributions. Here empty circles correspond to empty 
sites, while circles with small dots indicate particle positions. 
Parabolic curve at a site depicts full distortion of the lattice 
ions at that site with corresponding energy —g^ujo {+g^Loo) if 
the electron is present (absent) at that site. 
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FIG. 3. (Color online) n-particle normalized clustering prob- 
ability [NCP(n)] for number of sites =12, number of electrons 
= 6 (3 t,3 4-), t/ujo = 1.0, and when (a) g = 2.2 and (b) 
3 = 2.8. 



IV. RESULTS AND DISCUSSION 



The half-filled case (i.e., average concentration of one 
electron per site) is not very interesting because, owing 
to exclusion of double occupancy, we have one electron in 
every site and consequently no hopping occurs. Here the 
system behaves like a Heisenberg AF chain with NN re- 
pulsion having no effect on electronic correlations. How- 
ever, at quarter-filling there are enough vacancies for elec- 
trons so that, besides the spin interaction term, NN re- 
pulsion and different kinds of hopping in our model can 
contribute to a rich phase diagram. Therefore, we con- 
centrate on the ground state properties at quarter-filling 
and analyze the different phases exhibited by our effec- 
tive Hamiltonian in Eq. (|lip and the quantum phase 
transitions (QPT) between these phases for large values 
of g, i.e., in the strong e-ph coupling limit. 

We use a modified Lanczos algorithm^ to calculate the 
ground state. To get the basis states in the occupation 
number representation, we exclude double occupancy and 
permit only the remaining three possible electronic states 
on any site (i.e., t- i, and no particle). The size of the 
resulting Hilbert space prohibits us from calculating the 



ground state for large systems. For example, the number 
of basis states for 16 sites with 8 electrons (4'f, 44.) is 
^^Cs x*C4 = 900900 which is quite large but manageable 
for calculating the ground state. But for 20 sites with 10 
electrons (5t, 5J,), the calculation of ground state requires 
^°Cio C5 basis states which is prohibitive. 

We find evidence for several different types of phases 
depending on parameter values. Because there is a ten- 
dency to form phase separated clusters, various methods 
are used to identify the phases. To characterize phase 
separation, we use a "n— particle normalized clustering 
probability" parameter NCP(n) (defined precisely in Ap- 
pendix C) which actually measures the normalized prob- 
ability of rt— electron clusters in the ground state of the 
system. However, NCP(n) does not tell the number of t 
and 4 electrons or their arrangement in the cluster. To 
determine the different phases, we need to employ addi- 
tional measures of ordering such as spin-spin correlation 
functions, structure factor, etc. 
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FIG. 4. (Color online) Plot of the correlation functions (a) 
^"'{1) portraying an AF cluster for U/t= 15.2, t/uja = 1.0, 
and g = 2.2; and (b) W'"'{2) versus U/t again for t/cjo = 1.0 
and g = 2.2. 



A. AF Cluster to Correlated Singlets Transition 

Our numerical calculations were done for t/ujQ = 1.0 
and t/ujQ = 0.5. Since the results are similar for the 
two values of t/too considered, here we report only the 
calculations for t/ujo = 1.0. By varying g and U/t, we 
change the relative importance of J and V terms and get 
the different phases of the system. 

It is clear from Fig. ^a.) that, for t/wo — 1.0 and 
g = 2.2, for values of U/t up to about 15.6 the NCP(6) 
« 1.0 which implies that all the particles in the system 
form a single cluster. This is because the effect of spin 
interaction dominates over the NN repulsion. But upon 
increasing U /t further, the system undergoes a first-order 
QPT with NCP(6) discontinuously dropping to almost a 
zero value and NCP(2) concomitantly jumping abruptly 
to a value close to 1. The value of NCP(2) remains close 
to 1 up to J7/t « 20.0. Interestingly, for n=l,3,4,5 par- 
ticles, NCP(n) remains almost zero up to U/t k, 23.0. 
For all values of g in the range 2 < 5 < 3, when U/t 
is increased, the system transits discontiimously from a 
single cluster to a phase comprising of pairs of particles. 



This is verified for g = 2.8 as well in Fig. I3J6). Moreover, 
we notice that the correlated pairs phase persists over a 
broader window oi U/t values when g is larger. In fact, 
for all the system sizes considered (i.e., number of sites 
N = 8, 12, and 16), we find that the systems manifest 
a first-order QPT from an AF cluster of size N/2 to a 
correlated singlet phase as U /t increases. 

Although Fig. |3] demonstrates the QPT clearly, it 
does not give information about the structure/order of 
the phases. We will confirm the AF order by analyzing 
the correlation between density fluctuations of spins (ct 
and cr' separated by a distance I) through the correla- 
tion function W"' {I) = ^ " with 

j 

N being the system size. Here (nj) is the filling factor 
which for quarter-filling is 0.25. Fig. HJ^a) shows the cor- 
relation functions W^"^' il) for U/t = 15.2, t/ujQ = 1.0, 
and g = 2.2, i.e., away from the transition. This clearly 
shows that in the cluster, the dominant arrangement of 
the spins is antiferromagnetic with correlation decreasing 
with distance I. This is expected because of the trans- 
verse spin fiuctuation term in the Heiscnberg interaction. 
Next, in Fig. ^h), W'^'^{2) is shown as U/t is varied. 
There is a sudden jump in W^^''~(2) indicating that the 
system breaks up into two-particle clusters from an AF 
cluster. 

Now, we need to determine the structure of the phase 
with NCP(2) « 1.0 in Fig. H The important question is: 
what is the nature of the pairs? The answer lies in the 
correlation function plot [l) shown in Fig. ([5]) for 

U /t ^ 18.0 (i.e., deep inside the phase). Here, we see that 
the value of ly^^(l) is approximately its minimum possi- 



ble value of —0.0625 which occurs when j, ^( 



-1) 



0.0. This means that each pair is made up of two op- 
posite spin electrons. Furthermore, VFl^-'-(l) = 0.0625 in 
Fig. ([5]) which is expected because Vl^''~-'-(l) should have 
a calculated value of 0.25 x 0.5 - 0.0625 = 0.0625. Next, 
we notice that, for I > 2, W^^il) and have the 

same values. This means that, given a a a pair, we get 
not only an 'W. pair at a distance I with a certain proba- 
bility but also a J,t P^^ir ^it the same distance and with the 
same probability. Thus a pair of opposite spin electrons 
located at sites j and j + 1 can either be a singlet or a 



triplet with S^^^^^i — 0. When acted upon by the operator 
the singlet state yields eigenvalue -1 



■^total 

sjsf 



while the S^^^^^ = triplet state gives -1-1. Thus to know 
the nature of the spin pairs in the system, we have calcu- 
lated the quantity ^{{Sj'S~_^^ + SjSj]^-^)) and obtained 

j 

a value very close to —3.0 for 6 electrons (3 t and 3 |) for 
the whole range of 16.0 ^ U/t < 22.5. Thus the phase 
is entirely made up of singlets. In Fig. [5l we also notice 
that W^^{1) and W^^{1) show a peak at Z = 4 and slightly 
lesser values for I — 3 and I ~ 5 which is indicative of a 
CDW. We call this phase a correlated singlet phase. A 
detailed analysis of this phase at various filling factors 
will be presented elsewhere^. 

The transition from an AF cluster to correlated sin- 
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FIG. 5. (Color online) Plot of correlation functions W"" (l) 
in the correlated singlet phase for g — 2.2, U/t = 18.0, and 
t/ujo = 1.0. 
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FIG. 7. (Color onhne) Correlation functions W"^ (l) at 
t/uo = 1.0, g = 2.2, and for (a) U/t = 35.0 depicting AF 
order in one sublattice; and for (b) U/t = 50.0 displaying 
paramagnetic phase in one sublattice. 
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FIG. 6. (Color online) The phase diagram of the Hubbard- 
Holstein model at quarter- filling in the t/ojo = 1.0 plane. 



glets can be explained by invoking Bethe ansatz re- 
sults. Suppose we have a Hamiltonian of the type 

H = JJ2 ^j-^j+i + Then, for the cluster 

i j 

regime, Bethe ansatz yields energy/site = —0.443 J-|- 2-^ 
where as for separated singlets in the correlated singlet 
phase the energy/site = — J(0.375) + y. Thus the clus- 
ter regime prevails when — 0.443J-I- < — J(0.375) -I- 
i.e., V < 0.136J. Now, if in Eq. dTT]) we neglect all the 
hopping terms (which is true for large g), based on the 
above analysis, we obtain the condition for existence of 
the cluster phase to be {2V—j) < 0.136 J or equivalently 
U/t < S.bMg^ujo/t. This condition gives a very good es- 
timation of the actual phase boundary in the phase dia- 
gram given in Fig. [51 



B. Transition from correlated singlets to spins in 
one sublattice 



Here we discuss how, for a fixed g and larger values of 
U/t, the correlated singlets break up into spins occupying 
only one sublattice . It is clear from Fig. EJa) that this 
is a second-order phase transition at smaller values of g 
(such as 5 = 2.2) because NCP(l) [NCP(2)] increases 
[decreases] continuously as U/t is raised. On the other 
hand, this transition tends towards first-order at larger 
values of (7 [say g = 2.8 as shown in Fig. Elb)]. This tran- 
sition can be explained as follows. As U /t is increased, 
J gets reduced and NN repulsion starts dominating over 
spin interaction leading to the break up of singlets and 
single spins separating out to occupy one sublattice only. 



1. Transitions at smaller g 

For smaller values of g (i.e., 2 < g < 2.6), the QPT 
is not so sharp. In fact, for g = 2.2 [as shown in Fig. 
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FIG. 8. (Color online) Plot of the correlation function W'^^{2) 
versus U jt. 



probability of getting one t spin at any site, in a 12- 
site system, is 0.25. Now, as the spins are residing in 
only one sublattice, the probabilities of finding the re- 
maining 3 \ and 2 f electrons, in any of the remaining 
5 sites of the same sublattice, are 3/5 and 2/5 respec- 
tively. Hence, W^^{1) = 0.25(3/5) - 0.0625 = 0.0875 and 
W'^'^(0 = 0.25(2/5) -0.0625 = 0.0375 for / = 2, 4, 6, 8, 10. 
The above predicted values of the correlation functions 
W^^(V) and W^^iJ) match well with the calculated val- 
ues of these functions for Ujt = 50.0, g = 2.2 as shown 
in Fig. [71[b). Another interesting point is that there is 
a first-order phase transition from AF order to param- 
agnetic order in one sublattice of the system as can be 
seen from the sudden jump (when U jt is varied) in the 
correlation function W^''"~'-(2) depicted in Fig. |51 




-0.08 I ' ' ' ' 1 

1 2 3 4 5 6 

1 

FIG. 9. (Color online) Correlation functions W" (V) depict- 
ing spin disorder in one sublattice for Ujt = 50.0, t/ujo = 1, 
and g = 2.8. 



|3lja)], the singlet pairs and separated single spins coexist 
in the range 23 < U/t < 30. But when U/t is increased 
further, almost all the pairs are broken and the electrons 
occupy one sublattice only. The correlation functions 
W^^H) and W^^{1), depicted in Fig. [Tja) for U/t = 35.0 
and g = 2.2, show that the dominant structure has AF 
order in one sublattice (i.e, the spins arrange themselves 
ast — 4- — t — 4-— ••■)• This can be surmised from the 
fact that the correlation functions VFl^-'-(Z) and W^^{1) 
have dominant peaks at alternate even values of /. This 
spin arrangement is preferred by the system because it 
can gain energy due to AF interaction (J) by virtual 
hopping to NN site and returning back. For very large 
values of U/t, the spin interaction strength J becomes 
negligible and the electrons do not shown any spin or- 
der but still occupy only one sublattice. For this spin 
disordered situation, the values of W'^^{1) and W'^'^{1) 
for the system under consideration [i.e., 12 sites with 6 
electrons (3 ti 3 i)] can be calculated as follows. The 



2. Transition for larger g values 

For g > 2.6, when U/t is increased, the correlated sin- 
glet phase transits to a paramagnetic state (in one sub- 
lattice) directly unlike the earlier case for smaller values 
of g such as g = 2.2 (see Fig. This transition is of 
first-order character as shown in FigI3l^6) where NCP(l) 
and NCP(2) show simultaneous sharp jumps in opposing 
directions. 



V. THE PHASE DIAGRAM 

The phase diagram of the Hubbard-Holstcin model at 
quarter-filling in a twelve-site system for the t/ujQ = 1.0 
plane is shown in Fig. [51 This is obtained by producing 
plots of NCP(n) versus U/t for various values of g. Two 
such plots are shown in Fig. ^ for values of 5 = 2.2 
and g = 2.8. Then, from the NCP plots, we find phase 
transition points in the g — U/t plane to obtain the phase 
diagram. Furthermore, using selective choice of parame- 
ter values, we also found that the qualitative features of 
the phase diagram remained the same even for the larger 
sixteen-site system. 

In the phase diagram shown in Fig. |6l the system 
transits, for all values of g in the range 2 < g < 3, 
from an AF cluster to a correlated singlet phase discon- 
tinuously when U/t is increased. For g < 2.6, further 
increase in U/t drives the system continuously from a 
correlated singlet phase to a single-sublattice AF phase. 
Next, at even higher values of U /t and again for g < 2.6, 
the system jumps discontinuously from an AF order to 
a paramagnetic order in a single sublattice. Contrast- 
ingly, for g > 2.6, the system transforms directly from 
a correlated singlet phase to a paramagnetic phase by a 
close-to-discontinuous jump. 
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VI. CONCLUSIONS 



Appendix A 



In this appendix, we will demonstrate the validity 
of Eq. Let us assume a Hamiltonian of the form 

H — Hq + Hi where the eigenstates of Hq are all 
separable and are of the form |n, m) ~ \n)ei <8) \'m)ph and 
Hi is the electron-phonon interaction perturbation term 
of the form given in Eq. ([8]) . 



By analyzing the probability of occurrence of different 
cluster sizes and by studying various correlation functions 
that result from our effective Hubbard- Holstein Hamilto- 
nian, we deduced the ground state phase diagram (Fig. 
in]) at quarter-filling and in the non-adiabatic regime. 
In particular wc find that strong electron-phonon cou- 
pling stabilizes a correlated singlet phase that is charge- 
density- wave like, which would be absent in the pure 
Hubbard model. Notice that this phase is also quite 
distinct from a Peierls-like ( "bond-order") wave, being 
driven by onsite correlations alone. Our analysis and 
results should be of relevance to systems such as ful- 
lerides, lower-dimensional organic conductors, magnetic 
oxides, interfaces between oxides, where there is evidence 
for both strong Coulomb correlation as well as electron- 
phonon coupling. In higher dimensions, the concurrent 
NN spin-spin and NN repulsion interactions, could lead 
to even richer physics such as competition between or co- 
existence of charge ordering and superconductivity, non- 
Fermi liquid behavior, etc. 

I 



After a canonical transformation, we get 

H = e'^He-^ 

=Ho + Hi + [Ho + Hi,S]+^[[Ho + Hi,S],S].iAl) 

In the ground state energy, we know that the first-order 
perturbation term is zero. To eliminate the first-order 
term in Hi, we let Hi + [Ho,S] = 0. Consequently, we 
obtain 



{ni,mi\S\n2,m2) 



{ni,mi\Hi\n2,m2) 

{-^rti .nil ^^2,1712^ 



(A2) 



We now assume that both the exponentially reduced 
electronic NN hopping energy te~^ and the Heisenberg 
spin interaction energy J are negligible compared to the 
phononic energy wq which is true for large values of g. 



This implies that (i^m. 



Then, Eq. ()A2p simplifies to the form 

ph{mi\Hi\Tn2)ph 



ph{mi\S\m2)jjh 



ph 



(A3) 



Now, using Eqs. (jAl[) and (|A3|) . wc obtain 



ph{mi\H\m2)ph ^ ph{mi\Ho\m2)ph + '^'^[ ph{mi\Hi\m)ph ph{m\S\m2)ph - ph{mi\S\m)ph ph{m\Hi\m2) ph] 



ii\Ha\m2)ph + l-'^ph{mi\Hi\m)ph ph{m\Hi\m2)ph 



j^ph _ j^ph pph _ pph 
^m2 ^fh J-^mi J-^fh 



(A4) 



Appendix B: Projection onto singly occupied subspace 

In this appendix we evaluate each term in Eq. (jlOp by projecting out double occupancy. Every fermionic operator 
Cjcr is subject to the replacementcjcr Cj(j(l — rijg) so as to incorporate the action of the single subspace projection 
operator Pg. In the first two terms of Eq. ([TU]), the action of Ps implies that spin index a = a' . Thus, these first two 
terms correspond to the process where an electron hops to a neighboring site and comes back. Then, the projected 
form of the first term is evaluated to be 

Ps (c]+i^Cj',p2ct^Cj_^i„) Ps = (l- n^+m)c]+i„Cja{l - n^a){l - njff)c]^Cj+i^(l - rij+i^) 

= (1 - nj+is)cj^^^Cj„cj^Cj+i„{l - njg). (Bl) 

Similarly, for the second term we get 

Psc]^Cj + i^P^cl^j^^Cj„Ps = {l^j+la)c]„Cj+i„c]^i„Cj^il^j^). (B2) 
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The above two processes are depicted in Figs. [TJb) and[lja) respectively and add up to 



(B3) 



Next, let us consider the third and fourth terms in Eq. (|T0|). These two terms depict the process of an electron 
hopping to its neighboring site and then to the next-neighboring site. The process described by the fourth term is 
shown in Fig. [Ijc) and the third term is hermitian conjugate to it. Obviously, both neighboring and next-neighboring 
sites should be empty and here too we should take a = a' . On using the projection operator P,, these two terms yield 



3 „ . p2„t „ p 



= (1 - nj + lff)(l - ?7.jg-)(l - nj_lg-) c]_|_^^(l - ??.jo-)Cj_lcr + H.C. . (B4) 

Finally, we will consider the last two terms in Eq. (jlOp . These two terms represent the consecutive hopping processes 
where an electron of spin a at site j hops to its neighboring site j + 1 {j — \) followed by another electron of spin a' 
at site J — 1 (j + 1) hopping to site j. This successive hopping process may be termed a' a pair hopping. Here, the 
spin indices a and a' can be same or different. 

Now, the fifth term can be decomposed into two terms as follows of which : 

= (1 - nj+is){l - nj^i„)c]^Cj+i„c]_^^Cjs + (1 - ?ij+ia)(l - njs){l - nj^is)c]_^^{-nj„)cj+i„ (B5) 

The sixth term is the hermitian conjugate of the fifth term. Thus, upon adding the last two terms in Eq. (jlOp. we 
arrive at the expression 



)(1 -"j-la) 



H.c. 



+ {l - nj+is){l - njs){l - rij^is) c]^^^{-nja)cj-ia + Yi.c. (B6) 

The second term in Eq. (|B6P corresponds to the process depicted in Fig. WL,d) and the hermitian conjugate part of 
the first term is indicated by Fig. [T]^e). 



Appendix C: Definition of n-particle normalized 
clustering probability [NCP(n)] 

We are dealing with a system of N sites containing 
an even number iV/2 electrons. The ground state has 
equal number of ^ and | spin electrons. We get the 
basis states {(pi) in occupation number representation by 
populating N sites with N/2 electrons using all possible 
combinations with the constraint that each site can have 
only 3 possibilities (i.e., ti ii and no particle) with 
double occupancy being excluded. The ground state 
is obtained as a linear combination of these basis states: 
= ^ ai4>i where are the probability amplitudes. 

i 

Then, we calculate the n— particle Normalized Clustering 
Probability, i.e., NCP(n) from the ground state using 
the following procedure: 



1. Initialize the clustering probability CP(n) = 0.0 for 
n = 1 to N/2. 



2. Consider a basis state (pi with a corresponding co- 
efficient ai. 



3. Find an empty site (say j) in the basis state (pi. 

4. Start searching sequentially from site j + 1 onwards 
for occupied sites with index larger than j. 

5. Count the number of occupied sites until another 
empty site, say fc(> j), is reached. Since the size n 
of the unbroken cluster of electrons between the two 
empty sites j and k is given hy n — k — j — 1, add 
of to CP(k — j — 1). Then, again start searching 
for electrons from site k onwards until the next 
empty site is reached and again obtain the next 
unbroken cluster size ni. Similar to the previous 
case, add of to CP(ni). Continue the searching 
process for the whole system, i.e., from site j to 
site {j -t- N) with site (j -f N) being equivalent to 
site j. 

6. Repeat steps 3 to 5 successively for all the basis 
states (pi with corresponding coefficients to get 
CP(n) where 1 n N/2. 

7. Finally, by normalization, calculate NCP 
for n— particle cluster using the expression 

N/2 

NCP(n) = CP(n)/ X] CP(n) for 1 n ^ N/2. 

n=l 
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